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We consider magneto-transport in a disordered two-dimensional electron gas in the presence of 
a periodic modulation in one direction. Existing quasiclassical and quantum approaches to this 
problem account for Weiss oscillations in the resistivity tensor at moderate magnetic fields, as well as 
a strong modulation-induced modification of the Shubnikov-de Haas oscillations at higher magnetic 
fields. They do not account, however, for the operation at even higher magnetic fields of the integer 
quantum Hall effect, for which quantum interference processes are responsible. We introduce a field 
theory approach, based on a nonlinear a model, which encompasses naturally both the quasiclassical 
and quantum mechanical approaches, as well as providing a consistent means of extending them 
to include quantum interference corrections. A perturbative renormalization-group analysis of the 
field theory shows how weak localization corrections to the conductivity tensor may be described 
by a modification of the usual one-parameter scaling, such as to accommodate the anisotropy of the 
bare conductivity tensor. We also show how the two-parameter scaling, conjectured as a model for 
the quantum Hall effect in unmodulated systems, may be generalized similarly for the modulated 
■ system. Within this model we illustrate the operation of the quantum Hall effect in modulated 

systems for parameters that are realistic for current experiments. 
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INTRODUCTION 



'-^J \ The problem of electron motion in a disordered conductor in a periodic potential and strong magnetic field displays 
£^ ■ a rich combination of phenomena of both classical and quantum origin. The complexity arises from the appearance 
of two independent types of periodicities, of the potential and of the cyclotron orbits, which may interplay in com- 

i plicated ways. Perhaps the most striking effect is known as Weiss oscillations^, whereby very large oscillations in the 

resistivity are induced by even a weak periodic potential at moderate magnetic fields. At different magnetic fields or 
temperatures, other types of contributions to the resistivity become significant. For example, at very small fields, a 
positive magnetoresistance results from a classical mechanism of channeling of orbit a 2 i 3 i 4 . At larger magnetic fields, 
the resistivity develops Shubnikov-de Haas oscillations, which originate from the onset of Landau quantization and 
may still be affected by the periodic potentiali* 5 * 6 * 7 ^. At even higher fields, the integer quantum Hall effect^ (IQHE) 

(~~} \ becomes operative due to the contribution of quantum interference processes. 

■ Much of the above phenomenology has been thoroughly analyzed theoretically, by a variety of techniques including 

£SJ ' quantum mechanical approaches 5,6 involving diagrammaticsiiS*ii or solution of a quantum Boltzmann equation 8 , 
and quasiclassical approaches^iiii 2 ^. The experimental realization of these systems is also well advancedi* 2 * 5 * 6 *^*^ 
with precise confirmation of the theoretical predictions already possible. So far however a theory for the influence of 
quantum interference processes in a periodic potential and high magnetic field has been missing. Such processes lead to 
weak localization corrections to the conductivity which become significant at low temperatures (see e.g. Refsi2 3 j24i2Si2 6 i 
for reviews). At high magnetic fields, they are also responsible for the operation of the IQHESiSii, whereby the Hall 
conductivity becomes quantized at low temperatures. At the same time, the observation of the IQHE in systems with 
q | a weak potential modulation is well within current experimental capability^*^. 

This paper aims to fill this gap by showing how the standard theory for weak localization and the IQHE may be 
generalized so as to incorporate the presence of a periodic potential. We employ a field-theory approach based on a 
nonlinear a niodeliiiSiiSiSJiSi which is well established in the study of mesoscopic disordered conductors. In deriving 
the appropriate field theory, we are able to show how previous theoretical calculations of the conductivity, within both 
quasiclassical^A^ and quantum mechanical approaches 7 , may be recovered in a natural way within the field theory 
formalism. This enables us to extend previous theoretical results in a consistent way so as to include the influence of 
quantum interference processes. 

The experiments of Weiss et o/i employed weakly modulated two-dimensional (2D) electron systems of high mobility 
with a well-known period, a ~ 300 nm, much less than the mean free path, £ ~ 10 /im. Such samples were engineered 
using a holographic modulation technique, based on the persistent photoconductivity effect in GaAs/Al x Gai_ x As 
heterostructures at low temperatures ~ 4.2 K. The Weiss oscillations appear in only one component of the resistivity 
tensor, p xx , when the modulation is in the x direction. Furthermore they appear at magnetic fields B such that the 
cyclotron radius R c — vf/uj c (where lo c = eB/mc is the cyclotron frequency, — e is the electron charge and Vp is the 
Fermi velocity) satisfies the commensurability condition 2R C = (n — l/4)a, for integer n. Hence the oscillations are 
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periodic in 1/B. In addition they are relatively stable with respect to temperature, suggesting a quasiclassical origin. 

At even higher magnetic fields, such that the cyclotron radius is much less than the period of modulation (R c <C a), 
the quasiclassical theory predicts that the p xx component shows a large, nonoscillatory increase proportional to 
B 2 , leading to a strong positive magnetoresistance. This result has been confirmed in experiments 14 for which the 
temperature was kept deliberately high so as to avoid the intervention of the IQHE. 

A limitation of the quasiclassical approach is that it fails to account for the renormalization of single-particle 
properties, such as the density of states and scattering lifetime, by the strong magnetic field according to Landau 
quantization. Even in unmodulated samples, quantization leads to oscillations in the density of states with respect to 
magnetic fields, and hence to Shubnikov-de Haas oscillations in the resistivity. In samples with the periodic potential, 
the Shubnikov-de Haas oscillations start to appear at higher magnetic fields than the Weiss oscillations (for R c an 
integer multiple of the Fermi wavelength <C a) . A quantum- mechanical approach that does allow for such quantization 
effects has been provided by Zhang and Gerhardts 7 - (see also Peeters and Vasilopoulos 8 ); it is a diagrammatic treatment 
that generalizes the approach of Ando and co-workers^ to modulated samples. The quantum-mechanical approach 
is then capable of describing both Weiss oscillations and the Shubnikov-de Haas oscillations (also affected by the 
modulation) at higher magnetic fields. 

In principle, the calculation of weak localization corrections in a strong magnetic field (unitary ensemble) is possible 
even in the presence of a periodic potential, by a generalization of the diagrammatic approach of Zhang and Gerhardts 7 , 
but the procedure would be complicated (although the simpler orthogonal case has been examined diagrammatically 
for a periodic magnetic field22i31). Instead, the calculation of high-order diagrams is more convenient in the field- 
theory formalism and, furthermore, with this method the possibility exists of calculating contributions of diagrams 
to all orders by the renormalization-group technique. We show how the field-theory takes the form of a nonlinear 
a model with a topological termiLiLi2i22i2i. The effective Lagrangian is slightly nonstandard since it contains an 
anisotropy in the coefficients, corresponding to the difference between the longitudinal conductivities in the x and y 
directions, due to the periodic potential. 

The effect of weak localization corrections to the conductivity for unmodulated samples is accounted for by a scaling 
(one-parameter scaling) of the conductivity with the system size. As a first step we derive the analog of the one- 
parameter scaling for the conduct ivity^iSii^iSSi in the modulated system by means of a perturbative RG analysis of 
the effective Lagrangian. We then turn to the study of the IQHE, implementing a generalization of a two-parameter 
scaling, which has been conjecturedSiiSiiS as a model for the IQHE in unmodulated systems. We examine how the 
resistivity tensor should be affected at low temperatures by the IQHE in modulated samples for parameters that are 
realizable in actual experiments. We see, for example, how the Hall conductivity becomes quantized under scaling at 
low temperatures, while in the regions between the plateaus the longitudinal conductivities develop peaks of differing 
heights according to the anisotropy in the x and y directions. 

Model. In the following, we employ the Hamiltonian for the disordered conductor in a magnetic field and periodic 
potential in two dimensions: 



H = H -V(r), 

1 / eA\ 2 

H " = 2m V ~) +U o cos (v x )i (!) 

Here A is the vector potential, so that V x A = (0,0,-B), where B is the perpendicular, uniform magnetic field. 
Also Uq cos{qx) / (— e) is the periodic potential which is taken to be weak (Uq <C ep, where cf is the Fermi energy) and 
with modulation period a — 2n/q. V is the disorder potential, which we assume to be ^-correlated in space, with the 
associated scattering time r. We assume that the mean free path £ — vft greatly exceeds the modulation period a. 

The plan for the remainder of the paper is as follows. In Sec. [H] we describe the field theory approach and 
derive the effective Lagrangian for the system. In Sec. IIIII we show how the field theory provides the scaling of 
the conductivity tensor under changes of length scale due to the contribution of quantum interference processes, and 
hence a description of the IQHE in these samples. Section llVl concludes with a summary and discussion. 



II. FIELD-THEORY APPROACH 



In this section we describe a field-theoretical approach to the description of a disordered conductor in the presence 
of a periodic potential and a strong magnetic field. The approach is based on a diffusive nonlinear a model, which 
will serve as a tool for the analysis of quantum interference effects in the remainder of this paper. In section III Al 
we discuss this model and the relation to previous semiclassicaliSi and quantum mechanical 7 approaches. Technical 
details of the derivation are presented in section III Bl 
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A. Description of field theory approach 

The field-theory apparatus that we employ is by now well established in the study of spectral and wave-function 
statistics of disordered conductors. It is based on a functional integral expression of electron Green functions 
in the presence of disorder, and in the diffusive regime takes the form of a nonlinear a model in terms of a Q 
matrisiiiSiiS. In two dimensions, the field-theory technique provides a convenient means of calculating weak localiza- 
tion corrections to the conductivity, a task that may become very cumbersome by the more conventional technique of 
diagrammatics22iS&. The field-theory may also provide a resummation of diagrams to all orders by a renormalization- 
group procedur o 17 i 18 i 19 i 28 . 

We show below that the effective Lagrangian for the disordered conductor in the presence of a periodic potential 
and strong magnetic field takes the form 

C[Q) = ^J drStv{a xx (\7 x Q) 2 +a yy (V y Q) 2 -a xy r 3 Q[V x Q,\/ y Q}} , (2) 

where of- is the classical conductivity tensor of the system, in units of e 2 /h. The 8x8 supermatrix field Q(r) satisfies 
the nonlinear constraint Q(r) 2 — 1. The supertrace operation is defined in Rcf. 19 . 

In the absence of a strong magnetic field or modulation, a renormalization-group analysis of the first two terms 
in Eq. J5J) leads to the well-known one-parameter scalingi 7 *i2ii24 2 ^ of the conductivity with system size due to weak 
localization. 

The final term in Eq. (J2J is known as a topological term and appears in the theory of the IQHE proposed by 
Pruisken and collaborators 20 21 . The influence of the extra term does not appear within perturbation theory, but 
becomes evident only through a nonperturbative analysis. Such an analysis has been conjecturedSis^L^ in the form of 
a two-parameter renormalization-group procedure. The two parameters are the longitudinal and Hall conductivities 
that follow coupled scaling equations with respect to changes of length scale. 

While the validity of the two-parameter scaling has been vigorously debated (see, e.g., Refs. 34,35 ), it has remained 
a valuable guide to experimental- 7 " and numerical^ data for a number of years. In the following, we take a pragmatic 
approach by not contesting the validity of the two-parameter scaling and its relation to the proposed field theory. 
Instead, we assume its validity and explore how it may be generalized to take account of the periodic potential. 

The main difference of the Lagrangian from the unmodulated case is the anisotropy in the diffusion coefficients 
for the x and y directions. The Lagrangian, therefore, contains three, rather than two parameters, which scale together 
under changes of length scale. A simple scale transformation, however, maps the Lagrangian to an isotropic version 
for which the usual two-parameter scaling may be applied. 

This Lagrangian applies in the diffusive limit, that is, to configurations of the Q(r) field that vary on scales much 
longer than the mean free path. This allows the calculation of the contribution of low-momentum relaxational modes 
to weak localization corrections to the conductivity. Although momentum relaxation is diffusive on large length scales, 
electron motion on the scale of the periodic potential is ballistic, since the modulation period is much less than the 
mean free path. To arrive at the above Lagrangian it is therefore necessary to integrate over degrees of freedom 
corresponding to length much smaller than the mean free path. 

One way to do so is to start from a description of ballistic electron motion on length scales much less than the mean 
free path. Such a description is provided by the "ballistic a model"^2ii£ which, as the name suggests, generalizes the 
diffusive a model to the ballistic regime. Starting from this model we show in Sec. Ill Bl how the contribution from 
short length scales may be integrated out. The result is the Lagrangian of a diffusive a model that describes the 
interaction of diffusion modes on scales longer than the mean free path. 

When derived in this way, the effective Lagrangian takes the form of Eq. , but with the er coefficients replaced 
by afj(h/e 2 ). Here af? correspond precisely to the components of the quasiclassical conductivity tensor as derived 
from the Boltzmann equation for the modulated systemic. One finds that, in the quasiclassical approach, only the xx 
component of the quasiclassical resistivity tensor is affected by the modulation, all other components being the same 
as in the unmodulated case, so that pf x = po = (2e 2 b'D)~ 1 and p% c y = —p^. = u c Tpo (here D = v 2 f t /d is the diffusion 
coefficient in dimension d). For details of the solution of the Boltzmann equation, we refer the reader to Beenakker 12 . 
Here we only remark that for moderate magnetic fields, weak enough that the cyclotron radius is much larger than 
the modulation period (R c 3> a), but strong enough that lu c t 3> 1, the quasiclassical result simplifies to 

— = 1 + ^-( — ) ^cr) 2 R c q C OS 2 (R c q-^), (i? c g»l). (3) 

Equation. 10 demonstrates the oscillatory dependence of the resistivity known as Weiss oscillations. For strong 
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FIG. 1: Resistivity components as a function of magnetic field B for two different temperatures using the approximation 
scheme of Peeters and Vasilopoulos-. Corrections to the Hall resistivity are small in this regime-*-. Parameters are a =200 nm, 
U =0.3 meV, n el =3.4xlO n cm" 2 , h/r =0.011 meV. 



magnetic fields, such that R c <C a, a very large oscillatory increase in p xx proportional to B 2 is predicted and has 
been observed experimentally in, e.g., Geim et al^. 

However, for high magnetic fields (lo c t 3> 1), the quasiclassical results are reliable only at sufficiently high temper- 
atures (fcsT 3> h/r), while at lower temperatures the effects of quantization on single-particle properties such as the 
density of states must be taken into account by a quantum-mechanical approach. These effects are not taken into 
account in the derivation via the ballistic a model, since here the magnetic field is treated only as a weak pertur- 
bation. Indeed, once temperatures are low enough for weak localization corrections to be significant (ksT ~ h/i"), 
then quantization is already well established (as fcgT -C Hlo c , since uj c t 3> 1). Hence this method of deriving the 
Lagrangian is unfortunately of no use in describing the effects of quantum interference processes at high magnetic 
fields. 

An improvement on this situation may be found by proceeding instead along a second route to derive the effective 
Lagrangian, following more closely the original lines of Pruisken 20 . He showed how to derive a a model in the presence 
of a strong magnetic field, including the effects of quantization. The Lagrangian is applicable at high values of the 
Landau-level index n, such that the fluctuations of the density of states may be neglected. We show below in Sec. Ill Bl 
how this method may be adapted to describe a disordered conductor in a strong magnetic field and a periodic potential. 
This route bypasses the intermediate step of a ballistic er model, instead computing more directly the final diffusive 
Lagrangian. 

While the resulting Lagrangian takes again the form displayed in Eq. J2J, now the er° coefficients correspond 
precisely to a^h/e 2 , where cr^ u are the components of the conductivity tensor calculated in the fully quantum- 
mechanical approach of Refi Since the effects of quantization are included this approach correctly describes certain 
features that are observed in experiment but are beyond the quasiclassical approach. Some typical results for the 
resisitivity components (p^) calculated within the quantum- mechanical approach are shown in Fig. ^ 

As well as the appearance of Weiss oscillations in at moderate magnetic fields (up to ~ 0.4 T in Fig. [TJ, we see 
that both p^ and show strong, in-phase Shubnikov-de Haas oscillations at higher magnetic fields. These latter 
oscillations reflect those of the density of states. At weaker fields, for which the Weiss oscillations in p^ x are visible, 
weaker, out-of-phase oscillations in p^y also appear. Again, the latter oscillations in p^y are beyond the quasiclassical 
approach. 

The quasiclassical result for the change in p^ c x (for oj c t 3> 1 and R c q 3> 1) may be reproduced by the quantum 
mechanical approach, as long as the temperature is high enough for the thermal broadening to be greater than the 
separation of the Landau levels, ksT ;§> hw c (see, e.g., RefsA&). As the temperature is lowered, however, the region 
of applicability of the quasiclassical formula shrinks to progressively lower magnetic fields: The sinusoidal Weiss 
oscillations start to interfere at low temperatures with the Shubnikov-de Haas oscillations. Indeed, at vanishingly low 
temperatures, the quasiclassical result becomes unjustified for all magnetic fields. 

The Lagrangian J5J), supplied with the values of the conductivities determined by the quantum- mechanical approach, 
may then be used for a reliable description of quantum interference processes at high magnetic fields, and hence the 
operation of the IQHE. 
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B. Derivation of the Field Theory 



In this section we describe in more detail the derivation of the effective Lagrangian (J5J). In the field-theory approach, 
it is necessary to average the functional integral expression for products of Green functions over the disorder configu- 
rations. To perform the disorder average, several possible techniques exist, including the use of either supersymmetry 
or the replica trick. Although we may use either of these two approaches, we choose here the supersymmetry ap- 
proach, since it is free of certain technical problems regarding nonperturbative calculations that appear for the replica 
approach. The starting point for both routes is a functional integral representation of the partition function, Z , in 
terms of " superfields" iff, *S>: 

Z = J D[V]P[V]D[S>,y}ex-p-£{y 7 y, V), (4) 

C = i J dr^{r)(E - H - i(5A)*(r). (5) 

Here the superfields '4/ contain eight components, corresponding to fermion/boson, retarded/advanced and time- 
reversal sectors 1 ^. Also A = diag(l, —1) in advanced-retarded space space, S is a positive infinitesimal and H is the 
Hamiltonian Q up to the substitution A — > T3A, where T3 =diag(l, —1) in time-reversal space. The short-range 
random potential V(r) is taken to be Gaussian distributed, according to 



P[F(r)] = ie X p|-^ J drV(r) 2 Y 

where Z provides the normalization. 

Ballistic a Model. The first route is to derive a ballistic a model from the Lagrangian in Eq. (JSJ. Andreev et 
al*& have suggested a derivation of the ballistic a model by means of an energy averaging procedure, leading to a 
Lagrangian of the same form as that originally proposed by Muzykantskii and Khmelnitski: 39 . Following the lines of 
Andreev et a/^S, one treats the disorder field V as a perturbation in Eq. (JHJ . Performing an energy averaging on the 
clean Hamiltonian leads to a term in the Lagrangian that is nonlocal and quartic in the 'f, VP fields. The quartic term 
is then decoupled by a Hubbard-Stratonovich field, Q(r l7 r 2 ), and the 'J, # fields are integrated out. 

The resulting Lagrangian is then simplified by subjecting it to a saddle point analysis, and performing a gradient 
expansion around the saddle point. Note that, since at this stage we are neglecting the effects of quantization on 
single-particle properties, we do not include the renormalization of this saddle point solution by the magnetic field, 
instead we treat the magnetic field as a perturbation. 

Following a Wigner transform, and after applying a semiclassical 44 approximation, the Q matrix becomes a function 
of position, r, and the direction of momentum, n, where n is a unit vector. We do not repeat these steps here but 
refer to Refsi^iiS for further details. The resulting Lagrangian (ballistic a model) corresponding to the Hamiltonian 
© is 

C[Q n ] = ^J drStrj— / ^Qn(r)Q„,(r) - 2 J ^Ar 3 T~ ^(n, r)T n } , (6) 

In this equation C° is the Louiville operator defined by 

£° = u(r)n.V + uj c dc/, — sin cfyu' (x)d < f> (7) 
where <f> is the polar angle of n, v is the density of states per spin direction at the Fermi level, and 

v(r) = vf ( 1 — -=— cos(gx) 

is the electron velocity at the Fermi surface. The supermatrix Q n (r) satisfies the nonlinear constraint Q n (r) 2 = 1, as 
a result of which it may be parametrized by 

Q„(r)=T n (r)AT n - 1 (r). 

Due to the presence of the magnetic field, the Q n field satisfies unitary symmetry, that is, [Qn, r 3] = 0. The Q n ( r ) 
field also satisfies the further standard symmetries Q n = Q n = KQ^K, where Q n = CQ^C T . The matrices C and 
K, as well as supertrace operation, are defined in Refi 1 ^. 
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The Lagrangian © is in its current form too complicated for our purposes: it describes fluctuations of <3 n (r) on 
length scales smaller than the mean free path (and longer than the Fermi wavelength) and for any dependence on n. 
At the same time, electron motion on length scales much longer than the mean free path is diffusive. The strategy 
therefore is to integrate out the modes corresponding to fluctuations on length scales smaller than the mean free path, 
to produce a Lagrangian that describes the interaction of diffusion modes on larger length scales. 

To perform the integration, we isolate the the fluctuations of lowest mass (termed massless), corresponding to 
matrices Q n (r) that are independent of n and vary more slowly than the mean free path, £. We then integrate out all 
other fluctuations that preserve Q n = 1. A similar procedure has been carried out by Wolfle and Bhatli^i to derive the 
diffusive Lagrangian for a disordered conductor with anisotropic masses in the x and y directions (see also Refs42ii^ 
for further examples). 

To perform integration over the weakly massive modes, we write 



Qn = U(r)Q° a (r)U(v). 

Here the U(r) matrix, obeying the symmetries UU = 1 and U = KU'K, represents the massless fluctuations, while 
Qn( r ) contains the weakly massive fluctuations. The matrix Q n (r) may, in turn, be parametrized by its generator 
P n (r), for example, bjti^ 

Vn \l-%Pn 

Here P n is off-diagonal in retarded- advanced space and satisfies the symmetry P n = — P n . Since the fluctuations 
represented by the Q n matrix are weakly massive, it is sufficient to treat them within a Gaussian approximation: 
Q n = A(l + 2iP n — 2P% . . .). Inserting into the Lagrangian |JBJ|, we find 

irvh f , „ f f dn , „ „ n . . „ f f dn dn' ,„ . 
C - — I drSt *{-J ^Ar 3 P n £V,n)P n + -y --{I* - P n P n .) 

-2i J ^Ar 3 P n «(r)n.^| 

where $ = UWU, # = f/2[$,A]A. The next step is to integrate over matrices P n . This amounts to a set of 
Gaussian integrals with a linear term in P n , and hence may be performed by means of a shift of P n . 
The resulting Lagrangian has the form 

C\Q\ = 1^2 / drStr{aT x (V x Q) 2 +^ c y (V y Q) 2 - <r^r 3 Q[V x Q,V y Q}} , (8) 

It is written in terms of a matrix Q(r) = U(r)AU(r) that depends only on position and varies on length scales much 
longer than the mean free path (and hence the modulation period, a). The Q(r) field also satisfies the nonlinear 
constraint Q(r) 2 = 1, as well as the symmetries Q = Q = KQ^K. We see that the Lagrangian takes the usual form 
of a diffusive a model with a topological termSS, the only non-standard feature being the anisotropy of a%% ^ er^j • 
The coefficients <r^ c are defined by 



< ■2e^Qdr'^^v(r)v(r')n i n' j T(n,n',r,r')y CM 



where the averaging in Eq. @ is over one period of the modulation a (this averaging arises naturally in the Lagrangian 
(j2J since the $ matrix varies slowly on this scale). Also, F(n, n',r,r') is defined by the Boltzmann-like equation 

(£°(r, n) - C\ F(r, r', n, n') = 8{v - v')5{n - ri). (10) 

Here, the collision operator is given as 

C{Fm = ~(F(<l>)-J ff^)- 
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As we will show now the afj coefficients in the Lagrangian JHJ are the components of the quasiclassical conductivity 
tensor as derived by Beenakker— . In his approach, the Boltzmann equation is expressed^ (see also RefsiSi&iiii) in 
terms of the distribution function F(r,n) 

(c° - (?) F(r, n) = -e«(r)E • n, (11) 

and the current J(r) is given by 

T{v) = -2ei/(F(r,n)w(r)n 4 ), (12) 

where the averaging is over both the velocity direction n and one period of the modulation, a. The relation J 1 = o~fjEj 
defines the quasiclassical conductivity tensor cr^ c , whose inversion gives the quasiclassical resistivity tensor pfj. 
Comparing Eqs. ((10(1 and 1)11(1 . we see that r(n, n', r, r') is related to the distribution function F(r, n) by 



f dn' 

F(r,n)=-e / dr' u(r')E ■ nT(r, r', n, n'). 

J 27T 

Comparing also Eqs. JJjJ and i|12|) . we see that the <7^° coefficients in the Lagrangian (|SJ) calculated by this method 
coincide precisely with the components of the quasiclassical conductivity as derived 12 from the Boltzmann equation 
(fill) . In this way we see that we have redcrived the quasiclassical results within the field-theory formalism. 

As discussed above, the inadequacy of the Lagrangian JSJ is that it neglects the renormalization of the single-particle 
properties by the strong magnetic field. At high magnetic fields (oj c t » 1), once the temperature is low enough for 
the weak localization corrections to be significant (fc^T ~ h/r), the quasiclassical results for the conductivity^ have 
already become unreliable due to their neglect of quantization. In order to include these effects, and hence account 
for weak localization effects reliably, we need to follow a different route to derive the Lagrangian. Such a route for 
unmodulated systems has been provided by Pruisken2£, whose method may be adapted (as we show here) to include 
a periodic potential. This method includes the renormalization of the saddle point equation for the Q matrix by the 
strong magnetic field. It does not require a derivation of a ballistic u-model as an intermediate step, but provides 
more directly the final form of the Lagrangian. While the first route contains certain parallels with the quasiclassical 
approach 1 -, the second route contains closer parallels with the quantum-mechanical approach 7 . 

Generalization of Pruisken derivation. We now present the second route to the derivation of the effective Lagrangian 
(|2*|). starting again from form JSJ for the Lagrangian. The approach now is to average over the short-range disorder. 
This produces a term in the Lagrangian that is quartic in the 'J, ^> fields. The quartic term may then be decoupled 
by a (now local) Hubbard-Stratonovich field Q(r). After integrating out the ^> fields, one finds 



C[Q] 



«/r<!^StrQ 2 -istrln 
It 2 



-i -E + Uq cos(qx)) +SA + 



2t 



(13) 



where 7r M = ~ifV\7 M — er^A^/c. Here Q satisfies the symmetry Q — Q. 

The Lagrangian 1)13(1 . in principle, provides an exact description of the system, although in its current form it is 
too general to be useful. Instead, one proceeds by finding the saddle point value of Q that minimizes the Lagrangian, 
and by performing a gradient expansion about this minimum. The saddle point equation may be written as 



Q = — < r| (ir^ + U cos(qx) -E + i5A+ ^) |r > . ( I 4 1 



In order to find the matrix inverse of the operator in Eq. ((14(1 . we make use of the eigenvalues and eigenfunctions 
of the operator H — + Uo cos(qx). Using the Landau gauge A = (0, Bx, 0), one may write the eigenfunctions 

— 1/2 

of Ho in the form ^jjk, n {x,y) = L y exp(iT3ky)(f> nyXo (x), where L v is a normalization length. The center coordinate 
xo = l 2 B k remains a good quantum number despite modulation, where lg — (hc/eB) 1 / 2 is the magnetic length. The 
0n,x o ( a; ) are eigenfunctions of the Hamiltonian 

H 2 d 2 1 

H x0 = -^-^a + 2 mw c( x ~ x o) 2 + U cos(qx). (15) 
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In the absence of modulation (Uq — 0), Eq. I|15|) represents the Hamiltonian of a harmonic oscillator. We proceed 
by a first-order perturbative expansion modulation Uq. This turns out to be a very good approximation for typical 
parameters that we consider, such as a weak periodic potential. It is very difficult to improve on this approximation 
analytically, although exact numerical diagonalizations have been performed for the density of states^. Within the 
perturbative expansion, the eigenvalues e n (xo) — e n (xo + a) of the Hamiltonian (|I5|I are given by 

£n(^o) - E n (x ) = E n + U„ COs(qX ), (16) 

where E n — hcu c (n + 1/2) are the unperturbed Landau energies, and u n — Uq exp (-1/2X) L n (X), where X = q 2 l 2 B /2 
and L n are the Laguerre polynomials^. 

We see that the modulation lifts the degeneracy of the Landau levels and discrete levels are broadened into bands, 
whose width depends on the band index n in an oscillatory manner (due to the behavior of the Laguerre polynomials 
at large n). It is this oscillatory dependence that leads to the Weiss oscillations in the resistivity. The validity of 
the perturbation theory depends only on the smallness of the u n parameter, which is assured for large values of the 
Landau-level index n. We refer the reader to, e.g., Refsi£*£i& for further details of this perturbative expansion. 

We also make the assumption that in this basis, the saddle point value of Q is independent of the Landau indices 
and xo: this approximation is analogous to the C-number approximation (CNA) introduced by Zhang and Gerhardts 7 
for the self-energy matrix in the presence of modulation, and is valid as long as the magnetic field is not too high. In 
order to go beyond the CNA, one would need to generalize the saddle point Q to include a matrix structure in the 
space of Landau indices and a dependence on xo : such a task is of interest for future work but beyond the scope of 
this paper. 

The saddle point equation (|14fl may now be written 



= r^> - / dx (17) 

2r t"io E-E n (x )-i6A - V ' 



where Lq is the width of the Landau level in the absence of modulation and 



— fiuv — . 
2tt t 



(18) 



The saddle point equation i|17|) coincides with the self-energy equation derived in the self-consistent Born ap- 
proximation (SCBA) by Zhang and Gerhardtsi, under the replacement ihQ 1 ' 2 /(2r) — > T, R ' A , where indices refer to 
advanced/retarded space and £ denotes the self-energy. The solution for Q is then of the form iQ = eo + iApo, 
where po is proportional to the density of states. In the absence of modulation, the density of states reduces to that 
determined by Ando and co-workers^ within the SCBA (as shown by Pruisken2£). 

Having identified the saddle point value Q, we proceed by performing a gradient expansion of the Lagrangian i|13|) 
around the saddle point. We follow very closely the calculation of Pruisken^ii, although we use the supersymmetric 
rather than the replica formulation. We use the representation 

Q(r)=T(r)P(r)T- 1 (r), 

where the P fields are diagonal in retarded/advanced space and represent the massive modes. The procedure is to 
integrate over P and T separately. To do so, we split the fields of integration, 



DQ = J DP J £>Texp(Strln[/[P]]), 
in the process acquiring the associated Jacobian I[P] (as discussed in Ref ~&). £ becomes 

C[P,T] = ydr|-Strln[/[P]] + -^StrP 2 



--Strln 
2 



ihP 

i[E- T^tt^T - U cos(qx) - — 



To integrate over the P modes, wsplit C into two parts: 

C [P] = j-Strln[/[P]] + ^StrP 2 -istrln 



ihP 

E - 7r„7r M - U cos(qx) - — — 

At 



5C[P.T] = C[P,T]-C [P}. 
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Integration over the P fields then proceeds by cumulant-expanding 5C[P, T] with respect to £o[-P]. In turn, 5C[P, T] 
is computed by a gradient expansion of the Str In term up to second order in the combination = T _1 V ^T. The 
propagators for the latter expansion are of the form 

/ ihP V 1 

g(r,r') =< r| I E - ir M ir* - U cos(qx) - -^r( r )j k' >, 

weighted with respect to Co [P] . A typical second-order term is 

A yP = JdrJ dr'Str[^(r)][^(r')](g(r',r)7r^(r,r')^), 

where the averaging is over P with respect to Cq. At this point, we now exploit the assumption that the T(r) 
matrices vary in space more slowly than the modulation. Thus the averages of the products of the Green functions 
are short ranged with respect to the T matrices, and are translationally invariant after averaging over one cycle of 
the modulation. This allows us to perform a gradient expansion in the propagator averages: 



< 5 (r',r)7r^(r,rV) = K(r-r') 

= K^S{r ~ r') + K<&5(t - r')V^ u + ..., 

here the averaging includes that over a cycle of the modulation. A series of Ward identities may now be used as 
in Ref^S to simplify the resulting expressions up to second order in D^. The final Lagrangian is then of the form 
displayed in Eq. ©, where Q = TAT- 1 . 

Again, we have the standard form of a nonlinear er model with a topological term, although with the anisotropy 
a xx 7^ a yy Again, Q satisfies Q 2 = 1 and the symmetries Q — Q — KQ^K. Also the bare coefficients, er^ , are the 
components of the conductivity tensor at zero temperature in units of e 2 /h, according to the formulas 

-- f dr'(^,,( 5 1 (r,r')-5 2 (r,r')) 7 r x ,,( 5 1 (r',r)- 3 2 (r',r))), 

TO J 

0,1 i OJI 
xy "i xy ' 

J / dr'{{^ a '(r,r'^ t g 2 (r',v)) - (n, a Hr,r')w x!) Hr' »)} , (19) 

TO J 

here the averaging is over P as well as over one cycle of the modulation with respect to r. The propagator averages 
may now be evaluated by their saddle point values with respect to P. 

The bare coefficients tr^- are the components of the conductivity tensor, in units of e 2 /h, calculated in the SCBA at 
zero temperature. In the absence of modulation, the expressions (|19|) for the bare conductivity then reduce to those 
computed in the SCBA by Ando and co-workers 2 ^. In the presence of modulation, they reduce to the conductivity 
tensor computed in the SCBA by Zhang and GerhardtA of - = cr^ u (the equivalence of may be seen by following the 
working of Streda^l). Therefore, in contrast to those derived via the ballistic a model, the bare conductivities in the 
Lagrangian J2J represent the values calculated by a fully quantum-mechanical approach. 



7° 

xx,yy 



"xy 

-OJI 
"xy 



III. WEAK LOCALIZATION AND QUANTUM HALL EFFECT 

Having derived the final form of the Lagrangian (J2J , we are now in a position to use it to calculate weak localization 
corrections to the conductivity tensor. The effect of the corrections may be expressed as a scaling of the conductivity 
tensor under changes of length scale. A perturbative treatment of the Lagrangian J2Jl is valid in the limit of large 
longitudinal conductivity ((7® x ^> !)■ Within perturbation theory, however, the Hall conductivity is not renormalized 
by weak localization corrections. Instead, the effect of the topological term in the Lagrangian may only be made 
apparent through a nonperturbative analysis. Such an analysis has been conjectured 21 ' 31-32 in the form of a two- 
parameter renormalization-group procedure. The two parameters are the longitudinal and Hall conductivities which 
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follow coupled scaling equations with respect to changes of length scale. The derivation of the scaling equations 
presented in RefsiSiiS is based on approximating configurations of the field-theory parameter as a gas of instantons. 
While the validity of this derivation is being vigorously debated (see, e.g., Refsi ? 4 '?^ ), it has remained a valuable guide 
to experimental^ and numerical data for a number of years. In the following, we take a pragmatic approach by not 
contesting the validity of the two-parameter scaling and its relation to the proposed field theory. Instead, we assume 
its validity and explore how it may be generalized to take account of the periodic potential. 

A. Weak localization 

First we subject the Lagrangian to a perturbative RG analysis. The perturbation theory is valid in the limit of 
a xx 3> 1. The topological term in the Lagrangian does not contribute at any perturbative order and hence a xy is 
unrenormalized. For an unmodulated system, a xx — a yy and the perturbative RG procedure has been explained in 
detail, e.g., in Refp 1 ^. 

In the presence of modulation, the a model is slightly nonstandard due to the anisotropy of the longitudinal 
conductivities, a xx ^ a yy . However, quantum interference processes in a diffusive, anisotropic conductor have been 
considered before, for example by Wolfle and Bhatt^i, where the origin of the anisotropy was envisaged as due to a 
difference in effective electron masses in the two directions. In the latter paper, a diffusive a model was derived with 
the same form as Eq. although without the topological term. 

When a xx ^ a yy , the two parameters, a xx and <J yy , follow coupled flow equations under changes of length scale. 
While it is straightforward to derive and solve the two flow equations, it is also instructive to follow a different strategy 
whereby we perform a scale transformation after which the Lagrangian (J2J maps to an isotropic form: we scale 

= ^K>1) 1/4 , y' = y(^/<4) 1/4 , (20) 

under which the Lagrangian J5J) transforms to an isotropic a model: 

£ [3] = i i dr'Str{a°(VQ) 2 - a° xy T 3 Q[V x ,Q,V y >Q]} , (21) 



where a = y a xx a yy . The perturbative scaling for a may now be derived in the standard wayiS: the flow equation 
is (to leading order) 

da _ 1 
d\ogL 2ir 2 a 

By integrating the flow equation from microscopic to macroscopic length scales (of the order of the system size L), 
we find the conductivity 

1 /2 

^)=*°(l-^logf) . (22) 

Using the fact that the ratios a xxyy /a are invariant under the scaling, we may recover a xx , yy (L) from a(L) as 
follows: 

a xx (L) = a(L)^Ja° x /a° y , 

a yy (L) = a(L)^/a° y /a° xx . (23) 

We remark that the perturbative scaling derived from a field theory as above in the presence of anistropy in the 
conductivities in the x and y directions (due to an anisotropy in the electron mass) has been confirmed by direct 
diagrammatics (in the orthogonal ensemble) by Wolfle and Bhatli 4 ^. 

Equation H22|) shows how weak localization corrections affect the longitudinal conductivity within perturbation 
theory. We can see how the conventional scaling applies to a in the presence of the periodic potential. By inputting 
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FIG. 2: Plot of bare conductivities, providing a starting point for scaling, for a modulated vs homogeneous system. Parameters: 
a =140 nm, U =0.2 meV, T =0.048 PI meV and n ei =lxlO n cm 2 . 



the SCBA values for a xx and a yy weak localization corrections are obtained, according to Eq. they depend on 
both bare values a xx and a yy . Corrections to a xx , for example, are influenced by the strong oscillatory behavior of a yy . 
Upon inverting the conductivity tensor to find the resistivity tensor, the weak localization corrections will therefore 
give rise to additional oscillatory corrections to the resistivity tensor as a function of the magnetic field. As mentioned 
earlier, strong Weiss oscillations in p xx are accompanied by weak out-of-phase oscillations in p yy as long as quantum 
interference effects may be neglected. It is interesting to notice that weak localization corrections to p yy oscillate in 
phase with the (dominant) oscillations in p xx . 

In practice^!, rather than changing the system size L, one varies the temperature to change the effective system size, 
which is given by the diffusion length, (HD/^bT)) 1 / 2 . The perturbative result (|22|l is only valid when the corrections 
are much smaller than the bare conductivities, and hence may be difficult to verify experimentally. A potentially 
more promising strategy is to illustrate the effect of weak localization corrections when their contribution is relatively 
large such as is the case for the IQHE. 



B. Quantum Hall Effect 



We now generalize from the perturbative analysis to discuss the nonpcrturbative scaling that would follow from the 
Lagrangian (J2J. Again, we are able to make the scale transformation (|2(J[1 as above, which maps the Lagrangian to 
the isotropic version, Eq. 1|21[) . The coupled flow of the two coupling constants, a and a xy , may written in a general 
form, 



d^L = ^ a ^ jSiL-M*'*^- (24) 

The above beta functions, p and p xy , would take precisely the same form as the beta- functions that describe the flow 
of a xx and <j xy in the equivalent, unmodulated system. 

The starting point of the flow is determined by the bare conductivity tensor ay. In general, the coupled flow 
equations (|24l) for a and a xy need to be integrated up to length scales of the system size. The values of a xx (L) 
and o yy {L) may then be recovered from a{L) by Eqs. As a consequence of Eq. we see that the ratio 

o xx {L) / Oyy{L) remains constant under the scaling. This conclusion is furthermore independent of the form of the 
beta- functions, (3 and (3 xy , in Eq. J23J. 

In the absence of modulation, the SCBA values^ of the bare conductivities correspond to an approximate semicir- 
cular dependence of a® x on a xy . In the presence of modulation, the dependance of a on a xy is modified from the 
semicircle in a complicated way; a typical dependence, calculated using the scheme proposed in Refi, is shown in 
Fig.0 

The corresponding dependency of the bare conductivity tensor on magnetic field for the same parameters is part of 
Fig. below (thick curves). We remark that for the peaks centered around B=\.h T one may already expect vertex 
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FIG. 3: Evolution of the conductivity tensor for the modulated system with system size L under the two-parameter scaling. 
The thick curves are the bare conductivities (L = Lq), and the thin curves are the scaled conductivities for log(L/io) = 2. 
Parameters as in Fig. |5] 



corrections to become effective, although they are not included in the CNA approach that we have employed. Such 
corrections have the general tendency 1 ^ 1 - to enhance a yy in comparison to a xx . 

In order to provide an illustration of the IQHE in the modulated sample, it is necessary to assume a particular 
form of the scaling equations (|24|) . The following scaling equations2iiSi2i^& were derived (originally in the replica 
formulation) within a dilute gas approximation of the instantonic configurations of the Q matrix from the Lagrangian 
©: 

Tj— t = ~7rW - a 3 D c cos{2ira xy ) exp(-27ro-), 

c7 3 D c sm(2na xy ) exp(—2na). (25) 



rflogi 



The dimensionless constant D c is of order unity and is related to the density of instantons. It may be seen immediately 
from the form of Eq. (|25|l . that along the lines a xy = (n + 1/2), where n is an integer, the Hall conductance is 
unrenormalized. It may also been seen that the points (cr xy ,a) — (n,0), for integer n, are (attractive) fixed points. 
Upon scaling the system from microscopic to macroscopic length scales, the coordinates (<j xyi a) scale from the bare 
nonuniversal values towards the quantized values (n,0), for integer n. This tendency reflects the quantization of the 
Hall conductivity under scaling at low temperatures. 

As mentioned before, the validity of the approximations underlying Eq. (|25|l has long been the subject of debat o 34 ' 35 . 
Keeping this in mind we use these equations to provide an illustration of the operation of the IQHE in the modulated 
system in Figs. |3 and 0] 

Fig. |21 shows the evolution of the conductivity for the modulated system under the two-parameter scaling. We see 
how the Hall conductivity becomes quantized under scaling in the modulated system. Between the plateaus in the Hall 
conductivity, the longitudinal conductivities develop peaks under scaling of differing heights, due to the anisotropy 
in the system. Fig. ^ shows in more detail how these peaks in the longitudinal conductivities develop, while the ratio 
between the two conductivities remains constant under scaling. 

We remark that the behavior described above should be within current experimental capabilities: for example, 
Geim et aim* have studied the high magnetic field (R c <C a) regime, although they kept the temperature high to avoid 
the effect of quantum interference processes. Bagwell et alJ^> have also studied the IQHE in modulated samples; in 
their work however they have not focused on the regime of a weak periodic potential or effects that are independent 
of the device boundaries. 



IV. SUMMARY AND DISCUSSION 



In this paper we have considered transport properties of a disordered conductor in a periodic potential and strong 
magnetic field. We focused on the contribution of quantum interference processes, whose influence at high fields is 
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B [T] B [T] 

FIG. 4: Focus on the evolution of the longitudinal conductivities with system size L under the two-parameter scaling. Param- 
eters as in Fig- El an d log(L/Lo) = 0-4, from the broader to the narrower peaks. 



missing from previous approaches despite being responsible for a whole class of phenomena. For example, they lead to 
weak localization corrections to the conductivity and the operation of the IQHE at high magnetic fields. To this end, 
we introduced a field-theory approach, which is well established in the study of unmodulated disordered conductors. 

The effective Lagrangian of the field theory takes the form of a nonlinear a model, which describes the interaction 
of diffusion modes on large length scales. The presence of the strong magnetic field leads to an extra, topological term 
in the Lagrangian. The form of the Lagrangian is the same as for unmodulated systems, except for an anisotropy in 
the coefficients corresponding to the bare longitudinal conductivities in the x and y directions. 

We provided two different routes to deriving the Lagrangian. The first route was via a so-called "ballistic" a 
model, and demonstrated how the results of the quasiclassical approach^ may be recovered within the field-theory 
formalism. The drawback of this route, in common with the quasiclassical approximation, is that it neglects the 
renormalization of single-particle properies by Landau quantization. Consequently the resulting Lagrangian was too 
approximate to be useful in determining weak localization properties at low temperatures. The second route improved 
on this situation by including the effects of Landau quantization. It bypasses a derivation of a ballistic a model and 
instead follows more closely the original derivation of Pruisken^il for unmodulated systems. Indeed, while the first 
route contained more parallels with the quasiclassical approach, the second route contained more parallels with the 
quantum-mechanical approach » for the bare conductivity tensor. 

Having derived the effective Lagrangian, we showed how it leads to the scaling of the conductivity tensor under 
changes of length scale (and hence temperature). A perturbative renormalization-group analysis of the free energy 
leads to a generalization, to modulated systems, of one-parameter scaling for the longitudinal conductivities. Pertur- 
batively, the Hall conductivity is unrenormalized. In the regime of Weiss oscillations^, weak localization corrections 
give rise to an additional oscillatory dependence of the longitudinal conductivities as a function of the magnetic field. 
Due to their smallncss, these corrections may, however, be hard to detect experimentally. 

In order to describe the IQHE, whereby the Hall conductivity becomes quantized at low temperatures, a nonper- 
turbative analysis of the Lagrangian is necessary. This is provided by the conjecture of the two-parameter scaling 
for unmodulated systems. Assuming the validity of this conjecture and the underlying instanton gas approximation, 
we have shown how the two-parameter scaling may be generalized to the case of the modulated system. This has 
allowed us to illustrate the evolution of the resistivity tensor under the IQHE for parameters that are realistic for 
experiments. We find, for example, that the ratio of the two longitudinal components of the conductivity remains con- 
stant under scaling (an observation that does not itself depend on the assumptions used to derive the flow equations) . 
While the Hall conductivity still becomes quantized under scaling in the modulated system, between the plateaus the 
longitudinal conductivities develop peaks of differing heights due to the anisotropy in the system. 

There are several directions in which our analysis may be generalized. A simple generalization is to consider 
a periodic magnetic field, rather than potential, a situation which leads to a similar phenomenologj«^^&. The 
derivation of the field theory for this case follows very similar lines to the case of a periodic potential, with similar 
results. 

Another direction in which the analysis may be extended straightforwardly is to the study of the low magnetic field 
regime, in which a positive magnetoresistance has been observed^ and explained within a quasiclassical approach^**. 
As long as uj c t <C 1, one may derive the field theory for this regime according to the first route of Sec.|Hj neglecting 
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quantization; the coefficients of the field theory then coincide with the components of the conductivity tensor derived 
within the Boltzmann approach. The weak localization to the conductivity (or related quantities such as the dephasing 
time^Ji) may then be calculated using the Lagrangian for either the orthogonal or the unitary ensemble, or in the 
crossover between the pure symmetry classes. 

A less straightforward generalization is to improve on the CNA approximation of the quantum- mechanical approach 
and its analog in the field-theory formalism. According to this approximation, the self-energy matrix (and the saddle 
point solution of the Q matrix) are assumed to have a trivial structure in the space of Landau indices and coordinate 
xq. While this simplifies the analysis considerably, we have seen already how this approximation breaks down for very 
high magnetic fields, such that the cyclotron radius is much less than the period of modulation, R c <§; a. Improvement 
on the CNA is necessary, not only for such high fields, but also to study models of disorder with long-range correlations, 
for which vertex corrections are not negligible. Such models have been analyzed in the quasiclassical approach^ and 
shown to represent experimental data more closely in certain respects. 

Very recently it has been shown how to improve on the CN A 10 ' 11 within the quantum-mechanical approach, so 
as to analyze models at high magnetic fields or with long-range disorder correlations. In the field theory formalism, 
improvement on the CNA approximation would require enlarging the space of the Q matrix even further to include 
the additional matrix structure in the space of Landau indices, as well as a dependence on xq. Inclusion of such a 
structure within the Q field is to our knowledge a novel direction to pursue, although this task is left as a future 
project. 

A further area to explore is the case of a strong periodic potential, for which the first-order perturbative expansion 
in the potential that we use is no longer valid. While an analytical approach would be very difficult, a numerical 
treatment would be better suited to this regime. For a strong enough potential that U is of the order of u c , the 
smearing of the Landau levels is so great that they can no longer be individually resolved even at low temperatures. 
This leads to a quenching of the Shubnikov-de Haas oscillations (as noted by Beton et aim-) and hence one would 
expect the quantum Hall effect to be destroyed. 
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